Repulsion and attraction in high Tc superconductors 
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The influence of repulsion and attraction in high-Tc superconductors to the gap functions is 
studied. A systematic method is proposed to compute the gap functions using the irreducible 
representations of the point group. It is found that a pure s-wave superconductivity exists only 
at very low temperatures, and attractive potentials on the near shells significantly expand the gap 
functions and increase significantly the critical temperature of superconductivity. A strong on- 
site repulsion drives the A\g gap into a Big gap. It is expected that superconductivity with the 
A\g symmetry reaches a high critical temperature due to the cooperation of the on-site and the 
next-nearest neighbor attractions. 
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I. INTRODUCTION 



A few decades of study on high-Tc superconductor con- 
firmed the -y2-w&ve symmetry of the gap function at 
least in cuprate superconductors, relative to the conven- 
tional s-vifave gap symmetry [l], [3] • The newly discovered 
iron-based superconductors have a different gap sym- 
metry from cuprates, expressed as ^ cos /cj; cos fcy [1, |3] , 
named as s^-wave. Gap symmetries are believed to be 
results of point group symmetries of the superconduc- 
tor crystals. Tsuei and Kirtley provided a systematic 
description to the gap symmetries from the point group 
theory The relation between the gap symmetries and 
the interaction, however, are not clear in literatures. 
That is to say, what interaction determines a given gap 
symmetry is still unknown. 

The interaction in high-Tc superconductors is much 
more complicated than that in conventional s-wave su- 
perconductors. It is widely believed that the on-site 
repulsion of strongly correlated electrons drives Cooper 
pairs from s-wave to d-wave coupling. Monthous et al&\ 
gave a schematic picture for the interaction in cuprate 
superconductors as shown in Fig.l, where an electron 
locates at the center and another one hopping on the lat- 
tice. Due to the spin fluctuation repulsive potentials and 
attractive ones appear alternatively on the lattice sites 
but decay with increasing distances. A question is how 
they affect the gap functions of different symmetries and 
the critical temperatures of superconductivity. Under- 
standing this mechanism may help humans design new 
superconductors with higher critical temperatures. 

In this paper we first reveal the relation between the 
gap symmetries of a square lattice and the interactions 
on the lattice. This will provide a simple method for 
computing the gap functions of various symmetries. Us- 
ing this method we analyze the influence of the repulsive 
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FIG. 1: Potentials between a quasi-particle at the center of 
a square lattice (green spot) and another one moving on the 
lattice (red or blue spots). The radii of the spots represent 
the strength of the potentials^. Each of the four light green 
circles connects the neighbors on a shell. From inner to outer 
the shells are labeled as I, II, II, and IV. 



and attractive potentials on neighboring sites to the gap 
functions and critical temperatures of superconductivity. 
Finally we conclude that superconductivity with the Aig 
or Big gap symmetry may have much higher critical tem- 
perature than other symmetries, and an on-site attrac- 
tion favors a Aig gap symmetry but an on-site repulsion 
favors Big. 



II. GAP FUNCTION 

The BCS-type pairing Hamiltonian of superconductiv- 
ity is written as 

Hpair = T7 XI " ^)A\C^-kiC-k'iCk'\ (1) 
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where c\^,Ckc^,a =tji are the creation and destruction 
operators of electrons with spin up and down. Note that 
a factor has been added in the Hamiltonian. It is 
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necessary for convergence over k-space summation but 
missing in literatures for decades. The potential V^(q) is 
expanded into the real space 



nq)-(l/2)^nm) 
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(2) 



where V^(ni) is the potential energy between an electron 
at the center and another one at site m, as shown in 
Fig.l. 

The central site has the first shell neighbors, i.e., the 
nearest neighbors(n.n.), the second shell neighbors, i.e., 
the next nearest neighbors(n.n.n.), the third shell ones, 
and the fourth shell ones, etc. In general, only the on- 
site potential and the potentials on the four near shells, 
Vl, L — 0,1, 2, 3, 4, are important. Therefore, these five 
parameters determine the pairing Hamiltonian thus the 
properties of superconductivity. The gap function is 
given by the following well-known gap equation 7,] 



Ak = 



N 



where S,k = ^/e^ + is the energy of a quasi-particle, 
and ek is the free-electron energy on a square lattice 
given by et = — 2t(cos + cos ky) — At' cos k^ cos ky — ji. 
Changing the above gap function from k-space into 
the real space using a Fourier transformation Ak — 
^jjj A(m)e~*'' "^ one obtains 

A(m) = -IM y tanh(ek72fcr) 



(4) 



Compared to the case of k-space the above equation has 
much smaller dimensions, which is in general smaller 
than 21 x 21 as defined by the 20 neighbors of the four 
shells around the central one in Fig.l. The dimensions 
can be further reduced up to 5 x 5 when lattice sym- 
metries are included. In particular, obviously, a unique 
on-site potential Vq will result in a fc-independent poten- 
tial T^(k) = Vo thus the conventional s-wave constant 
gap function A(0). Dominant neighboring potentials 
Vi,V2,V^,Vi lead to the unconventional superconductiv- 
ity, such as the d-wave ones. 



III. THE SYMMETRIES OF GAP FUNCTIONS 
ON A SQUARE LATTICE 

The gap equation Q is symmetric about the point 
group G of the lattice since CRk = ^k,V{Rm) = 
V{m),R e G, that is 
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(5) 



This shows that the gaps A(m) transform according to 
the irreducible representations (IR's) of the point group 
G, i.e.. 



A(j)(i?m) = ^7^'J(i?-i)AW(m) 



(6) 



where T"{R^^) is the matrix of IR a for group element 
R~^. It is seen that gaps A(m) on the same shell dif- 
fer from each other at most by a matrix r"(_R^^). In 
particular, they differ only by a phase factor for one di- 
mensional IR's. This property is represented by the basis 
functions of IR's, i.e., gaps A(m) transform under the 
point group just as the basis functions do. 

A 2D square lattice has a point group symmetry of 
D4. This group has five IR's Aig, A2g, Big, B2g and 
Eg and their corresponding basis functions as listed in 
Table I. 



(3) Table I. Irreducible representations and basis functions. 



IR Basis function Allowed shells Wave type 
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+ y')/2 


0,I-IV 


s -1- s± 






IV 
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Big 


ix^~ 




I, III, IV 




B2g 


xy 




II,IV 


dxy 


Eg 


^x+iy\ 
\x — iy) 




I-IV 


p 



The first four of these IR's are one dimensional 
and the last one is two dimensional. The basis functions 
of the five IR's listed in Table I vanish except on the 
allowed shells. For example, a g-wave can only occur on 
shell IV, and a d^jy-wave on shells II and IV, etc. Then 



gaps A(m) are given by 



A(m) = ^^(^'"'^"') 

\jIR(Xm,y7n)\ 



(7) 



where fiR{xm,ym) is the value of the basis function at 
site {x„i,yjn), and A^^ is the magnitude of A(m) on 
shell Lm- 

Using ([7]) and making an inverse Fourier transforma- 
tion to A(m) one obtains the gap functions of every IR, 



A 



= Ao 
-4A2 



) -T- 2Ai (cos kx + cos ky) 
cos kx cos ky + 2A3(cos 2kx + cos 2ky) 
+ A A4(cos2kx cos ky + cos kx cos 2ky) 
Ayia, (k) = 4A4(sin k^ sin 2ky — sin 2kx sin ky) 
ABij,(k) = 2Ai(cosA:j: - cos ky) 
+ 2A3(cos2A:a; — cos2fcy) 
+ 4A4(cos 2kx cos ky — cos k^ cos 2ky) 
Asag (k) = -4A2 sin k^ sin ky 

— 4A4(sin 2kx sin ky + sin kx sin 2ky) 



(8) 
(9) 



(10) 
(11) 
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(k) = — 2Ai (iu sin — v sin ky) 

— 4 A2 (zM sin kx cos ky — v cos fcj; sin /cj^ ) 

— 4A3(iMsin2fca; — ■i;sin2A:j,) 

— 4A4[TO(sinfca; cos 2ky + 2 sin 2kx cos fcj,) 

— w(cos 2kx sin ky + 2 cos fc^^ sin2fcj,)] (12) 

where Al,L — 0,1,2,3,4 are the magnitudes of the 
on-site gap and the gaps on the four shells, and u — 
= "^(-i)- ^ s-wave gap Aq has been added 
into the gap function of Aig since the basis function of 
this IR is circularly symmetric. That is to say, a s-wave 
gap coexists only with that of the Aig symmetry. Which 
gap symmetry is selected by a superconductor is deter- 
mined by the lowest condensation energy of the system. 

These gap functions are simpler and clearer than the 
earlier results obtained by Tsuei and KirtleyQ. They 
gave only the first terms of these gap functions. We will 
see in the next section that the first terms may not be 
dominant. It should be emphasized that in general gap 
functions of different symmetries do not mix together, 
unless the symmetry ([5]) is broken. Mixing gap func- 
tions were considered due to minor anisotropy in the a 
and b directions on Cu02 planes, such as (1^2 _y2 -f idxy 
proposed by Ghosh et al^ and s + d^^-y^ by Miiller et 

IV. CRITICAL TEMPERATURES 

The gap functions of different symmetries are easily 
worked out through self-consistent computations to (jH) 
given potential energies Vl on different shells of a square 
lattice. It is found that the critical temperatures of super- 
conductivity of different gap symmetries are significantly 
different. 

In Fig. 2, a pure s-wave gap is created by a pure on-site 
attractive potential Vq = —0.5 and gaps with the Aig 
symmetry are created by the same on-site potential and 
other four —0.1, —0.2, —0.1, —0.01 on the four near shells. 
It is seen that the s-wave gap exists only below a very low 
temperature of about O.OOOlSi. Small attractive poten- 
tials on the near shells, especially on shell II, extends the 
on-site gap a few times and increases the critical temper- 
ature more than five times. Therefore, superconductivity 
is significantly enhanced by a pairing on the n.n.n. with 
the Aig symmetry. This pairing has two main compo- 
nents Aq and A2. This leads to a mixing s + s"*" gap 
function A(k) « Ao -I- A2 cos fc^, cos fcy. It has different 
signs at the center and (tt, 0) of the BZ just as the case in 
iron pnictides[12]. Further computations show that the 
potential on shell I, Vi, even if strongly repulsive, e.g. 
0.4, does not influence A(k) significantly. This is similar 
to the result obtained by Wang and co-workers [I^ for 
iron pnictides. This form of gap function is much more 
robust than the widely assumed form cos kx + cos ky in 
literatures [l3|. 

When the attractive potential on shell II, V2, increases 



5 3 



























b 

\ 

s-wave 


A, 

^2 







0.2 0.4 0.6 0.8 1 



x10-^ 

FIG. 2; Comparison between the gap of a pure s-wave and 
those of the Aig symmetry. The pure s-wave gap is created by 
a pure, on-site, attractive potential Vo ~ —0.5, and the gaps of 
Aig by the five potentials Vl = -0.5, -0.1, -0.2, -0.1, -0.01. 
Band parameters are set to be t = 1.0, t' = —0.125,/^ = 
—0.527 for a hole concentration of 0.125. 
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FIG. 3: Gaps with on-site potentials —0.1 and —0.5. The four- 
shell potentials are —0.1, —0.5, —0.2, —0.01. Band parameters 
are t = 1.0, t' = -0.125, ^ = -0.527. 



to —0.5 the gaps dramatically expand as shown in Fig. 3. 
Under the cooperation of attractive Vq and V2 the gaps 
increase significantly leading to a dramatic increase of Tc 
(about 52 times compared to the pure s-wave case). If the 
on-site attraction is weak the gap is mainly of a wave. 
Thus the superconductivity is significantly enhanced by 
the shell II pairing with the Aig symmetry. This pro- 
vides a possibility for researchers to increase the Tc of 
superconductivity by means of increasing the attraction 
on shell II atoms. This is just the case of the iron-based 
superconductors. It is expected that these superconduc- 
tors may have much higher critical temperatures than 
other series of superconductors. 

In cuprates, however, there is a strong on-site Coulom- 
bic repulsion due to the strong correlation effect. Thus a 
on-site pairing is energetically expensive in these super- 
conductors. The gaps with the Ai and Bi symmetries 
under different on-site potentials are shown in Fig. 4. It is 
seen that the on-site repulsion strongly reduces the gaps 
with the Aig symmetry but holds those with the Big 
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enhanced by increasing attractive potentials on other 
shells. It is found through computations that the po- 
tential on shell II, no matter attractive or repulsive, does 
not affect the Big pairing, but V3 has a weak influence 
to the gap, as shown in Fig. 5. Ai increases only slightly 
with more attractive V3 but even a repulsive V3 does not 
break the Big pairing. 

Gaps with other symmetries A2g, B2g and Eg are ver- 
ified through computations to be small and result in low 
critical temperatures. Superconductivity with the Aig 
and Big symmetries is most expectable to go into the 
higher temperature region. 



FIG. 4: Gaps under different on-site potentials Vo- Potentials 
on the four shells are -0.3,-0.5,-0.2,-0.01. t = 1.0, t' = 
-0.125, /i = -0.527. 



CONCLUSION 
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FIG. 5: Gap Ai with the Big symmetry vs potential on shell 
III. Other potentials are -0.1,-0.5,-0.2,-0.01 except V3. 
t = 1.0, t' = -0.125, 11 = -0.527 



symmetry. Therefore, under a strong on-site repulsion a 
pairing with the Big symmetry, i.e. a d^2_y2 wave, is 
much more energetically favorable. This is just the case 
in cuprates which display a d-wave superconductivity. 
A question is whether this Big pairing can be further 



In this work we studied the influence of repulsion and 
attraction in high temperature superconductors to the 
gap functions and transition temperatures. We proposed 
a systematic method to compute the gap functions us- 
ing the irreducible representations of the crystalline point 
group. Then we analyze the gap functions with various 
symmetries for different potential energies at different 
temperatures. It is found that a pure s-wave supercon- 
ductivity exists only at very low temperatures, and at- 
tractive potentials on the near shells significantly expand 
the gap functions and increase dramatically the critical 
temperatures. Especially the Aig gap is expanded by an 
on-site attractive potential and that on the next near- 
est neighbors, thus increases the critical temperature for 
tens of times. A strong on-site repulsion, however, blows 
up the Aig gap but keeps a Big gap survive. It is ex- 
pected that superconductivity with the Aig symmetry 
has a high critical temperature. 
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